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Abstract 

We report the first observation of a narrow charmonium state produced in the exclusive decay 
process — > K ± tt + tt~ J/ip. This state, which decays into ir + ir~ J /if), has a mass of 3871.8 ± 
0.7(stat) ±0.4(syst) MeV, which is very near the Md + Md* mass threshold. The results are based 
on an analysis of 152M BB events collected at the T(4S) resonance in the Belle detector at the 
KEKB collider. 

PACS numbers: 14.40.Gx,12.39.Mk,13.20.He 
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INTRODUCTION 



A major experimental issue for the cc charmonium particle system is the missing (i.e., 
still unobserved) charmonium states that are expected to be below threshold for decays to 
open charm and, thus, are narrow. These include the n — 1 singlet P state, the h c (lP), 
and, possibly, the J PC = 2~+ r} c2 and J PC = 2 1 3 D 2 ^ 2 (1D), the n = 1 singlet 

and triplet spin-2 D states, which are narrow if their masses are below the DD* threshold. 
The observation of these states and the determination of their masses would provide useful 
information about the spin-dependence of the charmonium potential. 

B meson decays are a prolific source of charmonium particles. The large B meson samples 
produced at 5-factories provide excellent opportunities to search for the missing states. 
The Belle group recently reported the first observation of the r] c (2S) via its KgK^-n^ decay 
channel in exclusive B — > KKsK~-k + decays based on an analysis of 44. 8M BB events 
Strategies for finding the remaining missing states have been presented by Eichten, Lane 
and Quigg 0; they note that a narrow ^ 2 (1D) should have a substantial decay branching 
fraction for the 7T + 7T~J/ifj final state. In this paper, we report on a search for a narrow 
peak in the 7i + n^J/ip mass spectrum from exclusive B + — ► K + -k + -k~ J/tp decays 3] using a 
152M BB event sample. The data were collected in the Belle detector at the KEKB energy- 
asymmetric e + e~ collider, which operates at the T(4S) resonance (v/i = 10.58 GeV). KEKB 
is described in detail in ref. Q; Belle is described in detail in ref. p. 



5-MESON SIGNAL SELECTION 

We use events that satisfy the Belle-standard hadronic-event selection requirements and 
have a pair of well identified oppositely charged electrons or muons with an invariant mass 
in the range 3.076 < M £ + e - < 3.116 GeV, a loosely identified charged kaon and a pair of 
oppositely charged pions. In order to reject background from 7 conversion products and 
curling tracks, we require the 7i + n~ invariant mass to be greater than 0.4 GeV. To reduce 
the level of e + e~ — > qq (q = u, d, s or c-quarks) continuum events in the sample, we also 
require i? 2 < 0.4, where i? 2 is the normalized Fox- Wolfram moment 0, and | cos#b| < 0.8, 
where 9b is the polar angle of the S-meson direction in the center-of-mass system (cms). 

Candidate B + — > K + 7i + 7i^J/ip mesons are reconstructed using the energy difference 

AE = E c B ms — E^ m and the beam-energy constrained mass M bc = ^(-^b^m) 2 — G !" 8 ) 2 ) 
where E^ m is the beam energy in the cms, and E B ms and are the cms energy and 
momentum of the B candidate. The signal region is defined as 5.271 < M^ c < 5.289 GeV/c 2 
and \AE\ < 0.030 GeV, which corresponds to ±3cr and ±2.5cr from the central values, 
respectively. 



RESULTS 

Figure ^a) shows the AM = M(tt + tt~£ + £~) — M(£ + £~) mass distribution for events in 
the AE-M\y C signal region. Here a large peak corresponding to if)' — > -K + -K~J/ip is evident 
near 0.59 GeV. In addition there is a significant spike in the distribution near 0.77 GeV. 
Figure \Bb) shows the AM distribution for a large sample of generic BB Monte Carlo (MC) 
events. Except for the prominent ip' peak, the distribution is smooth and featureless. In 
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FIG. 1: Distribution of M(7r+7r~ •£+■£") - M{l + f) for events in the AE-M hc signal box for (a) 
Belle data and (b) generic BB MC events. 



the rest of this paper we use M(7r + 7r~ J/ip) determined from AM + Mj/^, where Mj/^ is 
PDG [7| value for the J/ip mass. 

In Fig. |2l we show the M\, c distributions for 5-MeV M(7r + 7r~ J/ip) bins for the AM region 
around the spike at 0.77 GeV. Here the \AE\ < 0.03 GeV requirement has been applied. 
Significant I?-meson signals appear in the M{tt + t[~ J/ip) = 3870 MeV mass bin and its 
immediately adjacent bins, but nowhere else. The curves in the figures indicate the results 
of fits to the distributions using a Gaussian function to represent the signal and the ARGUS 
function for the background [8J|. Here the peak and width of the signal Gaussian are fixed 
at values determined from fitting M^ c distributions in the if)' peak region. 
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FIG. 2: Mbc distributions for events in the AE signal box for 5-MeV M(tt + it~ J ftp) mass bins in 
the vicinity of the AM = 0.77 GeV spike. The mass values listed in each panel correspond to the 
center of the bin. The curves are the results of fits described in the text. 
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Fit results 
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FIG. 3: The data points show the B-meson signal yields from fits to the individual M]-> c distribu- 
tions in the (a) M(tt + i{~ J/ip) = 3870 MeV and (b) if)' regions. The curve represents the results 
of the fit described in the text. 



The S-meson signal yields from the fits to the individual Mt, c distributions are shown 
as the data points in Fig. 0(a). The curve is the result of a fit with a Gaussian signal 
function and a second-order polynomial to represent the background. The fit gives a signal 
yield of 34.4 ± 6.5 events, a peak mass of M% eas = 3871.4 ± 0.7 MeV, and a width of 

a = 3.4 ± 0.5 MeV. The signal significance is estimated to be ^—2 \n(C / C max ) = 8.6a, 
where C max and £ are the likelihood values for the fitted signal yield and for a fit with 
the amplitude of the signal function fixed at zero, respectively. For comparison, we show in 
Fig. E(b) the corresponding distribution and fit for the if>' peak region. Here the fit results 
are: 456.6 ± 22.4 events; Mft eas = 3685.6 ± 0.2 MeV; and a = 3.4 ± 0.2 MeV. The fitted 
widths for the if)' and signal peaks are in good agreement; the fitted if)' mass is 0.4 MeV 
below the PDG value of M?P G = 3686.00 ± 0.09 MeV. 



Mass determination and limits on the total width 

We determine the mass of the signal peak by referencing to the well measured if)' mass: 

M x = MY as - Mft eas + M^P G = 3871.8 ± 0.7 MeV. 

Here the error is statistical only. Since we reference Mx to the precisely known PDG value 
for the if}' mass, the systematic error is small. The My measurement, which references the 
J /if) mass that is 590 MeV away, is within 0.4 MeV of the nominal PDG value. Variation of 
the mass scale from My to the Mx requires an extrapolation of only 186 MeV and, thus, 
can safely be expected to be less than this amount. We conservatively assign 0.4 MeV as 
the systematic error on the mass. 

The fitted signal width is consistent with the resolution. To determine a 90% confidence- 
level (CL) upper limit on the total width, we repeat fits to the distribution in Fig. Efa) 
using a resolution-broadened Breit-Wigner (BW) function to represent the signal. From 
the variation in —2 ln(£) with the input BW width we determine a 90% CL upper limit of 
r < 3.5 MeV. 
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Di-pion mass distribution and product branching fraction 



The open histogram in Fig. Ufa) shows the tc + it~ invariant mass distribution for the 
M = 3872 MeV signal region; the shaded histogram shows the corresponding distribution for 
events in the non-signal AE-M^ region, normalized to the signal area. The 7r + 7r~ invariant 
masses tend to cluster near the kinematic boundary, which is around the p mass; the entries 
below the p are consistent with background. Figure Efb) shows the tt + tc~ mass distribution 
for the ip' events. This distribution also peaks near the upper kinematic boundary, which in 
this case is near 590 MeV. 
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FIG. 4: The M (tt + tt^) distributions for events in the (a) M(ir + -K^ J/ip) = 3870 MeV signal region 
and (b) the ip' region. The shaded histogram is the sideband data normalized to the signal-box 



area. 



Assuming equal acceptances for both channels, we determine a ratio of product branching 
fractions for B + -> K+X, X -> n+n-J/ip and B + -> K+tp', ip' -> tt+tt- J/if> to be 



B{B + -> K+X) x B(X -> TT+vr- J /if}) 

B(B+ -»• K+t/}') X B{lp' ->• 7T + 7T- Iff) 



0.075 ± 0.014(stat) ± 0.007(syst). 



Here the systematic error includes possible different efficiencies for the M(7r + 7r~) > 0.4 GeV 
selection requirement for the two channels, estimated by comparing the ip' yield with and 
without this requirement. 

DISCUSSION 



The mass of the observed state is higher than potential model expectations 
of-gravity (cog) of the l 3 D c j states: the Cornell and the Buchmiiller-Tye 
both give M cog (lD) = 3810 MeV, which is 60 MeV below the observed state, 
state of the 3 D c j triplet, one expects the 3 D C 2 to be near the cog.) Moreover, 
has been identified with the ^"(3770). Identification of the observed state 
state would imply a 3 D C 2 — 3 D c x splitting of ~ 100 MeV. On the other hand, 
due to effects of coupling to real DD and virtual DD* states may reduce 
potential model calculations for these states 12l . 
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The 7r + 7r~ invariant masses for the M = 3872 MeV signal region concentrate near the 
upper kinematic boundary (Fig. Ufa)), as is also the case for if)' — > n + n^J/ip (Fig.HJb)). For 
the M = 3872 MeV signal region, however, the boundary corresponds to the mass of the p. 
The decay 3 D c2 — > p J /if) is an isospin- violating process, and these are strongly suppressed 
in if)' — > J I if) transitions. With more data, we could determine whether or not the 7r + 7r~ 
system is coming from p meson decay, which could have important consequences for the 
interpretation of this state. 

The measured mass of the state is within errors of the D°D*° mass threshold (3871.3 ± 
0.5 MeV Q). This would be expected for a loosely bound DD* multiquark "molecular 
state," such as proposed by De Rujula, Georgi and Glashow in 1977 

CONCLUSION 

We have observed a strong signal (8.6<r) for a state that decays to -k + -k~ J /ip with 

M = 3871.8 ± 0.7 (stat) ± 0.4 (syst) MeV 
T < 3.5 MeV . 

This mass value is about 60 MeV higher than potential model predictions for a ID charmo- 
nium state and equal, within errors, to M D o + M D *o. This coincidence with the D°D*° mass 
threshold suggests that this may be a DD* multiquark state. 
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